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Abstract 

A method is proposed for arbitrarily engineering the high-order harmonic generation phase 
achieved by shaping a laser pulse and employing xuv light or x rays for ionization. This renders 
the production of bandwidth-limited attosecond pulses possible while avoiding the use of filters for 
chirp compensation. By adding the first 8 Fourier components to a sinusoidal field of lO^^W/cm^, 
the bandwidth-limited emission of 8 as is shown to be possible from a Li^^ gas. The scheme is 
extendable to the zs-scale. 
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I. INTRODUCTION 



High-order harmonic generation (HHG) is the key technology for attosecond science. The 
shortest pulse durations currently achieved are below 100 as [IHS]- Highly efficient sources 
are available having a bandwidth of hundreds of electronvolts [1] being large enough to 
produce pulses of only 10 as duration — in case they can be generated without chirp in the 
future. Many properties of HHG radiation emitted from a gas target can be understood 
by studying a single atom. The three-step model [5j is the simplest model to describe the 
single-atom dynamics. The process starts when a strong laser field liberates the electron of 
an atom and subsequently drives it in the continuum. If ionization happenend at the right 
time, the electron can be accelerated back towards the parent ion after the field has changed 
its sign and can recombine along with the emission of an energetic photon. 

One prominent feature of HHG is that the emitted light has an intrinsic chirp, the so- 
called attochirp [6l [7]. The origin of the attochirp can be understood from the classical 
electron trajectories in the laser field: trajectories with different energies recollide at different 
times. Due to the widespread classical recollision times in a usual sinusoidal laser field, the 
emitted harmonic pulses have a longer duration than their bandwidth limit, that is the 
minimum pulse duration for a given spectral bandwidth reached when the spectral phase 
of the harmonic pulse is constant. To compress the emitted pulse, dispersive elements 
of either chirped multilayer x-ray mirrors [8], thin metallic ffims [U HU] or gaseous media 
[21 [H] are employed and even the use of grating compressors is attempted JTT\ . However, 
these techniques suffer from losses, rely on the material properties and the chirp is not 
well-controlable in these schemes. Additionally, it is required to select radiation from the 
positively chirped short trajectory via phase-matching before the compensation element. 
Other approaches [UlIS] try to exploit the reduced time window of recollision when increasing 
the wavelength of the driving laser at a constant bandwidth and field maximum. However, 
this way, it cannot be taken advantage of the increased cutoffi Without selection of a certain 
trajectory, it is possible to partially reduce the attochirp as shown in [HI [15] by adding a 
weak second-harmonic or subharmonic field. The attochirp problem is expected to get more 
significant and demanding when photon energies spanning far into the soft x-ray domain are 
reached in the future. 

In this paper, we propose a way to engineer the attochirp in a determined manner by 
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altering the harmonic generation process. In terms of the wave function, the following 
scenario is realized: the wave function localized in the binding potential is continuously 
freed resulting in a large spread in space and momentum. The electronic quantum dynamics 
is tuned in a way that after a certain time of propagation, the wave function spatially re- 
compresses at least along the propagation direction of the wave packet. It has its minimum 
width exactly at the time of recoUision but with a large energy bandwidth gained during 
propagation in the continuum. We show that when the laser pulse is shaped by adding a 
small number of low-order harmonics and employing soft x rays for ionization, attosecond 
pulses with arbitrary chirp can be formed including the possibility of attochirp-free HHG 
and bandwidth-limited attosecond pulses. 

The benefit of assisting the HHG process in a strong laser field with a weak high-frequency 
field has been demonstrated mainly for the purpose of enhancing the single-atom yield [161 - 
[20] , improving phase matching |2Dj or suppressing the relativistic drift [2T| [22] . On the other 
hand, femtosecond pulse shaping has been used to shape the HHG spectrum [23], to increase 
the HHG cutoff or for relativistic HHG [251 [2S]- Here, we employ x rays to ionize the 
electron with non-zero velocity and combine it with femtosecond pulse shaping to control 
the spectral phase of the harmonic spectrum. 

II. CLASSICAL ANALYSIS 

It is known |27] that the HHG process can be analyzed by taking only a few quatum orbits 
into account that correspond to classical trajectories. Therefore, a classical analysis is able 
to give first insight into our idea. In the first part of our paper, we start out by considering 
classical trajectories in a tailored laser field [see Fig. [l] (a)] and find the condition when 
trajectories ionized at different times recoUide at the same time. For the considered laser 
intensities, the classical trajectories have only a component along the polarization direction 
of the laser field, i.e. 1-dimensional trajectories are shown. We describe the principle of our 
method by discussing two example trajectories marked by a and (3 in Fig. [l](b) being ionized 
separated by a small time span We first focus on the point in time when trajectory /3 just 
starts. At that time, a has already been driven slightly away from the origin. The distance 
between the trajectories is 6xi ~ PiSti where pi is an initial momentum e.g. mediated via one- 
photon-ionization (atomic units are used throughout unless stated elsewise). The momentum 
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FIG. 1: Schematic of the recolhsion scenario: a) A half cycle of the tailored laser field (black). The 
red line is the assisting x-ray pulse, b) Different classical trajectories in the field of (a) which start 
into the continuum at different times but revisit the ionic core at the same time. 



difference at that time is given by 6pi ~ —EidU because a has already been decelerated by 
the laser field Ei. From now on, the momentum difference 6p = 6pi is conserved during 
the whole propagation time. The separation of both trajectories at recollision after a time 
T is therefore given by 6xe ~ Sxi + 6pT = {pi — EiT)6ti. 6xe = reflects the condition 
on the initial momentum and on the electric fleld at ionization necessary for simultaneous 
recollision of the two trajectories: 

Ei ^ pi/r (1) 

The previous discussion applies only to the case of two single trajectories ionized with 
an inflnitely small time separation but is sufficient for demonstration of our principle idea 
of simultaneous recollision for trajectories ionized from a flnite time window. 

In Fig. [l| we chose the simplest possible laser fleld to fulflll condition ([T]): a flrst plateau of 
strength Ei and duration Atj which is followed by a second higher plateau of duration Atu. 
In this procedure, the fleld strength E^ of the second plateau is chosen flrst. It determines 
r for trajectories /3. Then, the flrst plateau is found having a small slope determined by op- 
timizing a polynomial expression to have simultaneous recollision of all classical trajectories 



4 



(approximately condition Q). 

Crucial for the previous scheme is the ionization with a non-zero velocity. The condition 
can be reached by using x rays of frequency ux which are responsible for ionization. In 
this case, the initial momentum is pi = \/2{uj^^-^Jp) where Ip is the binding potential. The 
X rays co-propagate with the laser field and are spectrally filtered out before the harmonic 
radiation reaches the detector, ujx can be chosen in a way (with the appropriate laser field 
according to Eq. ([T])) to meet a resonance of a suitable absorber at ujx- Alternatively, the 
X rays could propagate with a tiny angle to the propagation direction of the laser. In Fig. 
[!} we only consider trajectories ionized within At/ and having a direction pointing upwards 
to the potential (starting in positive direction in Fig. [T]). The initial momentum could also 
be directed downwards but in this case, the classical electron would not recoUide and no 
harmonics would be emitted. A second branch of trajectories that are ignored in the figure 
are those trajectories emerging from the second plateau during At//. They re-encounter the 
core region but all at different times and would break the desired scenario of simultaneous 
recoUision. Therefore, ionization during the second plateau has to be suppressed which 
demands for the use of an x-ray pulse [sketched as red wiggeld line in Fig. [T] (a)] rather 
than a cw x-ray field. The x-ray pulse must have a non-vanishing field strength only during 
At/ being of order of 1 fs. Moreover, as medium, we choose a gas of ions because their 
high ionization potential is necessary to suppress tunnel ionization. The plasma could be 
generated via laser ionization with a strong pre-pulse. 

We briefly comment on the energy distribution of the recoUiding trajectories. The classical 
trajectory marked by /3 has a distinct role because it experiences approximately a constant 
laser field. Therefore, it is symmetric to the turning point and also recollides with same 
momentum pi as it started. Trajectories starting prior to /3 recollide with a higher energy 
but at the same time as /3 assuming the first plateau-like structure is chosen in agreement 
with Eq. g as it is in Fig. [TJ Trajectory (3 plays the role of the trajectory with the 
lowest energy ux- The duration of the first plateau At/ determines the velocity difference 
Ap ^ EiAtj and the energy bandwidth Auq ~ ^(Ap + Pi)'^ + Ip — = \^'P^ + ^VVi- 

In the optimum case shown in Fig. [T| all trajectories recollide simultaneously. However, 
under real experimental conditions, deviations of the laser and x-ray field from the optimal 
conditions result in a non-zero time window Atg of recoUision. The classical model em- 
ployed a discrete x-ray frequency ux rather than a finite bandwidth as a real pulse. The 



5 



maximum allowed bandwidth can be deduced from the former model. The x rays ionize the 
electron with an initial velocity and arrange this way the initial displacement 5xi between 
two trajectories. If ujx deviates from its optimal value by ^ux-, the initial momentum will 
deviate by 5pi ~ 6u)x/{'2pi) and result in a additional displacement 6xi = 6pi6ti which is not 
compensated for. Therefore, the final wave packet size is of order 6piAtj resulting in a time 
spread of 

Atr = Sp.AtJp = (2) 

2p,p 

which has to be smaller than the envisaged pulse duration. 

III. STRONG-FIELD APPROXIMATION MODEL 

So far, purely classical dynamics was considered. In order to model the single-atom HHG 
yield, we use the strong-field approximation (SFA) and include the laser field within the 
dipole approximation. This way, the Fourier transformed dipole matrix element is given by 
[191I22] 

dt / dt' / d='q($o(t)|x|p + A(t)/c)Ejf(t')(p + A(t)/c|x|$o(t'))e"^^''^P'*'*'^ (3) 

oo J ~oo J 

where ^^{t) is the ground state wave function of the employed zero-range potential, 
5'q(p,t,t') = + A.{t") /cY /2 + Ip}dt" +uJxt' — q^t the classical action and q the harmonic 
number. In the long wavelength regime the highly oscillating integral can be evaluated by 
using the saddle-point approximation. Expression dsj) is approximated by a sum 



= E W T^('^o(te)|x|p. + A(t)/c)Ex(t.)(p. + A(t)/c|x|<|.o(t.))e"*^'(P-*-*') (4) 
^ Y <lei{Ss) 

over the saddle points s = {ps,te,ti) defined by 

'\ps + A{t")/c]dt" = (5) 



[p, + A{U)/c]y2 = ujx-lp (6) 
[p, + A(te)/c]V2 + Ip = qco (7) 

for S(ij) = didjS where i,j e {px,Py,Pz,'t,t'}. In case ux — Ip is positive all saddle points 
are real and, therefore, tj and are the ionization and recoUision times, respectively, being 
also a solution of the classical equations of motion. 



We analyze the condition required for (near-) bandwidth-limited emission of high har- 
monics. The time-dependend intensity of the emitted light bursts is given by 
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(8) 



where the phase Re Sg{ps, t^, ti) ^ + oiiiuOq — ujc) + \oi2iojq — + ... is crucial for the 
duration of the attosecond burst. We Taylor-expanded Re Sqijps^te^ti) about the central 
harmonic frequency Uc- As outlined in [6j,i7j the linear coefficient of the expansion, the group 
delay, is simply given by 

d 



du, 



-Re Sq{psj ^ej ti 



— Re teioJq 



(9) 



because the saddle-point equations lead to vanishing partial derivatives. The linear chirp, 
the group-delay dispersion (GDD), is given by 
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being in first approximation responsible for the duration of the harmonic pulse. Hence, 
the pulse can be considered to be bandwidth-limited if the quadratic term in the Taylor- 
expansion fulfills the following demand |^Q;2(Aa;g)^| ^ 2ti. This yields a criterion for the 
classical recollision time window Atg resulting from the bandwidth of the trajectories: 



Ate < 



An Ait 



Atr 



a 



(11) 



with the bandwidth-limited pulse duration Atp = a/Auq and the parameter a of order of 
unity determined by the spectrum. This refiects that only an approximate condition At^ ~ 
is necessary for bandwidth-limited HHG emission. In the previous section, we discussed in 
detail how the condition of simultaneous recollision of classical trajectories Atg ~ can be 
f unfilled. We will show that it is sufficient to add only a few low-order Fourier components 
of the fundamental laser frequency to the laser field in order to achieve an optimized field 



which fullfills Eq. (11). Moreover, Eq. (11) evidences that the chirp compensation becomes 



more difficult for larger bandwidth. 

In our classical and quantum descriptions, we used the common assumption [T6H20] of 
neglecting the x-ray field for the continuum propagation of the electron. This means drop- 
ping the term Sx{p,t,t') = /j!(ir([p^ + A{t)/c] Ax{t)/c + ^AI(t)/c^) in the definition of 



Sq{p,t,t'). We derive an approximate condition for the maximum electric field strength of 
the X rays being allowed that this approximation holds and that the HHG pulse duration is 
not influenced by the x-ray field. When evaluating at the previously determined saddle 
points, the additional shift of the recoUision times by the x-ray field is: 

dUg Oti OUJq Ote OUJq OPs OUq 

^ -[p + A{U)/c] A^{U)lcp- ^ -^1^ (12b) 

OOJq y/uj]^ OUq 



In Eq. (12a) we dropped the second term because ^ is negligible due to the simultaneous 
recoUision and the third term because it contains an integration over the heavily oscillating 
A^{t)/c. Thus, the overall additional recoUision time spread 

Ate,x = ^^<Ate (13) 
y/u]c AUq 



has to be much smaller than the recoUision time window to justify in our case the neglection 
of Sx- 

In summary, bandwidth-limited HHG emission can be reached when the recoUision time 



window determined by the classical dynamics in the laser field fulfills Eq. (11) and the x-ray 
field satisfies Eq. ^ and Eq. (13). 



IV. GENERATION OF ATTOSECOND PULSES 



In the following, we describe an implementation of the scheme. We start from an optimal 
field shape determined in the way described above with a fundamental frequency of u = 
0.06 a.u. Then, we represent the field as a Fourier series and only take the Np = 8 lowest 
frequency components of its spectrum into account. Following this procedure, the pulses in 
Fig. [2] a) (solid black line) were found. Its relevant classical trajectories ionized by a one- 
photon transition are shown as solid red lines. We can observe that te{ujq) is approximately 
constant and, therefore, we can expect bandwidth-limited harmonic emission. We compare 
it to the traditional case of a sinusoidal laser field with similar parameters as in [1] (dashed 
lines). Evaluating the dipole moment within the SFA [Eq. (|4])], we find the attosecond pulse 
generated by the fields in Fig. |2] (a). The results are shown in Fig. [2] (b). The duration of 
the pulse generated by our method (solid line) is 86 as at full-width half maximum (FWHM) 
only being slightly longer than a pulse generated by a sinusoidal field in combination with a 
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FIG. 2: a) Classical electron trajectories (red lines) are shown for two different field configurations 
(thick black lines): The dashed curves are for a conventional sinusoidal laser field whereas the full 
lines are for a shaped laser field consisting of Np = 8 harmonics that is assisted by continuous wave 
x-rays {dx = 60 eV, Ix = 9 x 10^^ W/cm^) triggering ionization, b) Attosecond pulses emitted 
from the laser fields proposed in a) . The solid line shows the intensity for the case of the proposed 
scheme whereas the dotted and dashed lines correspond to the attosecond pulses for a conventional 
field with a 40 eV bandwidth either with or without phase compensation, respectively. 



perfect chirp compensation (78 as) (dotted line). The dashed line shows the uncompensated 
pulse with a duration of 130 as. In all three cases, the same frequency window between 60 
eV and 110 eV was employed and the x-ray field {ux = 60 eV, J^; = 9 x 10^^ W/cm^) was 
chosen to have the same average ionization rates as by tunnel ionization in the sinusoidal 
field. The two pulses from the sinusoidal pulse were obtained from Neon using only the 
short trajectories whereas in our case He^ has to be employed. The laser frequency is in 
all cases u = 0.06 a.u. Although the single-atom emission rate of both examples are on the 
same order, the photon yield of our setup will be lower because the gas density is restricted 
by a maximum value of 5 x 10^^/cm~^ as determined later. We estimate an emitted photon 
number per half cycle of 10^ in the traditional case (density 10^^/cm~'^) and 10^ in our case 
from a volume of 200 yUm x 200 x 1 mm. 

So far, the attosecond pulse production with durations little below 100 as was considered. 
The time difference of the uncompensated (~ 130 as) and compensated (~ 80 as) pulses 
were moderate. Now, we discuss our proposal for future experiments where the bandwidth 



is of order of keV and Eq. (11) imposes a stronger constraint, demanding for a larger phase 
compensation over a larger frequency spectrum. In the following, we present two examples 
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with production of pulses below 10 as and 1 as, see Fig. [3j 
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FIG. 3: a) shows the laser field composed of only 8 Fourier components that illuminates an Li^"*" 
atom with Ip = 4.5 a.u.. The parameters of the additional x-ray field are indicated in the first row 
of Tab. |l} The created 8 as pulse is shown in b). c) displays the laser field needed to create a pulse 
of 800 zs duration from Be^"*" atoms with Ip = 8 a.u. and the parameters indicated in the second 
row of Tab. |l| The respective pulse is shown in d). The dashed red lines are the temporal phases 
of the pulses. 



The laser field strength and the x-ray frequency have to be increased compared to the 
example before to obtain a much larger bandwidth and are indicated in Tab. [!} This way, 
attosecond pulses with a FWHM of 8 as and 800 zs are formed. The pulses are almost 
bandwidth-limited as can be seen from the constant phase (red dashed line) in the main 
part of the pulse. Here, the linear phase term given by the central frequency of the pulse 
was subtracted from the phase. 

However, the photon yield emitted from the target is very low for several reasons: In 
the two examples, ux ^ Ip is valid because a much larger initial momentum is required, 
however, the ratio renders ionization inefficient. We limited the x-ray intensity according to 



Eq. (13). On the other hand, one-photon ionization exhibits a larger angular distribution 
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TABLE I: Parameters for the two examples in Fig. [Sj Np represents the number of Fourier 
components contained in the fundamental pulse, Ip its peak intensity, lox the x-ray frequency 
employed for ionization. Ex its field strength, Aug the achieved HHG bandwidth, Tp the HHG 
pulse duration, pmax the estimated maximum gas density and Nph an estimate of the emitted HHG 
photon number per half cycle emitted from a volume of 200 fim x 200 /xm x 1 mm having the 
maximum density. 

resulting in a large spread of the electron as soon as the initial momentum becomes larger. 
Finally, the gas density is limited to a small value which will be discussed shortly. 



V. MACROSCOPIC EFFECTS 



In the remainder of this paper, we discuss consequences of applying the scheme to a 
macroscopic gas target. Due to dispersion, the initially optimal pulse shape will be deformed 
during propagation. We estimate the impact of the dispersion by investigating the pulse 
shape after 1-dimensional propagation through a plasma of length L and refractive index 
n„ = w 1 — where pe is the electron density. In this case, each Fourier component of 



frequency Ug propagates with the phase velocity Vph = c/n and an analytic expression can 
be obtained: 

lit) oc I 5]]^2i+i 4-+ie-'"^^+^f*+("^~"^^+^^^/'l ' (14) 
j 

The influence of the different atomic transition lines of the medium on the dispersion is 
omitted because the free electron background forms the largest contribution to the disper- 
sion. When calculating the driving pulse shape of Fig. |3] a) after a propagation length of 
L = 1 mm, we find a maximum ion density of 5 x lO^^/cm^ and 7 x lO^'i/cm^ to maintain a 
duration of the harmonic burst below 10 as and 1 as, respectively. Similarly, we can specify 
the precision of the phase of the different Fourier components that is required. The allowed 
fluctuations of the different components in terms of time delay is of the order of 25 as and 
2.5 as, respectively, in agreement with the time delays caused by the plasma dispersion 
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discussed before. The sensitivity is lower for harmonics with lower energies. 

Apart from causing a deviation from the optimized pulse form, dispersion can also lead to 
phase mismatching. Due to the rather low gas density, we do not expect a dramatic phase 
mismatch. To achieve phase matching, we propose either to exploit the geometry of the 
laser focus or to use quasi-phasematching schemes as employing a weak counterpropagating 
IR field |28l EH], weak static fields [30] or modulated wave guides [3TH33] . 

So far, the recoUision scheme was discussed in the spotlight of bandwidth-limited emis- 
sion of attosecond pulses. The scheme can also be applied as a new type of pump-probe 
technique where the atom is excited or probed at a precise time by the recoUiding electron. 
The recollision time can be controlled via shaping the driving pulse. Moreover, the spec- 
tral diversity of the simultaneously recolliding trajectories is an excellent condition for the 
observation of continuum-continuum harmonics [3l] . 

We thank R. Moshammer and C. Buth for fruitful discussions. 
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